SYMMETRIES OF THE BLACK-SCHOLES EQUATION 



PAUL LESCOT 



Abstract. We determine the algebra of isovectors for the Black-Scholes equa- 
tion. As a consequence, we obtain some previously unknown transformations 
on the solutions. 

MSC 34 A 26, 91 B 28 



1. Introduction 

We study the algebro-geometrical structure of the Black-Scholes equation. 

The computation of the symmetry group of a partial differential equation has 
proven itself to be a very useful tool in Classical Mathematical Physics (see e.g. 
|Harrison~E stabro ok 1971] ^ . as well as in Euclidean Quantum Mechanics (see e.g. 
[Lescot-Zambrini 2004] and |Lescot-Zambrini 2008] h It was therefore natural to 
try and use the same method in Financial Mathematics. 

After setting the general framework (§2), and performing some preliminary re- 
ductions (§3), we determine (§4) the isovectors for the Black-Scholes equation in a 
way broadly similar to the one used for the backward heat equation with potential 
term in the second of the two aforementioned joint papers with J.-C. Zambrini. 
Our computation turns out to suggest Black and Scholes' original solution method 
( Black-Scholes 1973 ) of their equation; in particular, the quantities r — %- and 
r + ^j- appear naturally in this context. As corollaries, we determine (§5) the struc- 
ture of the Lie algebra of the symmetry group of the equation, then (§6) we obtain 
some interesting transformations on the solutions. 
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2. Generalities and notations 
We shall be concerned with the classical Black-Scholes equation : 

dc l 2q2 d 2 c dc 

for the price C (t, S) of a call option with maturity T and strike price K on an 
underlying asset satisfying S t = S (see [Black-Scholes 1973] , where a is denoted by 
v, C by w, and S by x). As is well-known ([Black-Sc holes 1973] . p. 646), the same 
equation is satisfied by the price of a put option. 
We assume a > 0, and define 

T 2 



and 

ft is useful to remark that 



a 

7 '~ r ~Y 



o 2 



r z s z 
2^2 +r = 2^2 



We intend to determine the isovectors for (£), using the method applied, in [Harriso n-Estab rook 1971] . 
pp. 657-658 (see also |Lescot-Z ambrin i~2004] . pp. 189-192) to the heat equation, 
and in [Lescot- Zambrini 2008] . §3, to the (backward) heat equation with a potential 
term. 

Let us set x = ln(5) and 

<p(t,x) :=C(t,e*) = C(t,S) ; 
then <p is defined on R + x R. One has 

dC_ _ ld<p 
~dS ~ Sdx~ ' 
d 2 C _ 1 dtp 1 d 2 ip 

as 2 " - ~52^ + s^dx^ 

and 

dC_ _ dip 
~dt ~~dl ' 

Equation (£ ) is therefore equivalent to the following equation in p : 
dip a 2 d 2 ip dip dip 
-d-t + Y { d^-d- X )+r d- X - riP = H£lh 

that is : 

dp a 2 d 2 ip „dip 
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3. Computation of the isovectors : preliminary reductions 

dip dip 

Let us set A — — — and B = — — , and consider thenceforth t, x. ip, A and B 
dx dt \ . . 

as independent variables . Then (£2) is equivalent to the vanishing, on the five- 
dimensional manifold M — R+ x R 4 of (t,x,ip, A, B), of the following system of 
differential forms : 

(3.1) a = dip - Adx - Bdt , 



(3.2) da = -dAdx - dBdt , 
and 

(3.3) /3 = (B + fA- rip)dxdt + -<T 2 dAdt . 

Let I denote the ideal of AT*(M) generated by a, da and f3 ; as 
d/3 = (dB + fdA - rdip)dxdt 

(3.4) = da(dx - fdt) + a(-rdxdt) e I , 

I is a differential ideal of AT*(M). By definition (see Harrison-Estabrook 1971]), 
an isovector for (£ 2 ) is a vector field 

(3.5) N = N t!L + N x® +N <plL + N AjL + N BjL 

at ox dip oA oB 

such that 

(3.6) C N {I) C I . 

Using the formal properties of the Lie derivative ([Harrison- Estabrook 1971] . p.654), 

one easily proves that the set Q of these isovectors constitutes a Lie algebra (for 

the usual bracket of vector fields). 

In order to determine G, we may use a trick first explained in 

[Harrison- Estabrook 1971] . p. 657, that applies in all situations in which there is 

only one 1-form among the given generators of the ideal / (see also [Lescot- Z ambrini~2 8 , 

p.211). 

Let N G Q ; as Cjy(I) C J, one has £jv(a) E / =< a, da, f3 >, whence there is a 
0-form (i.e. a function) A such that C^la) = Xa. Let us define 

(3.7) F:=N\a = N v -AN X -BN* . 
This can be rewritten as 

(3.8) Xa = C N (a) =N\da + d(N\a) = N\da + dF , 
whence 

(3.9) N\da = Xa-dF , 
that is 

(3.10) N\ (-dAdx - dBdt) = Xa-dF 
i.e. 

(3.11) -N A dx + N x dA-N B dt + N t dB 

= Xa-dF 

= X(dip - Adx - Bdt) - dF . 
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Whence (letters as lower indices indicating differentiation, as usual) 



(*) 



-N B 
—N A 

= A — F,„ 



-XB 
-XA 



F r . 



N x = -F A 
^ N* = -F B ■ 

Using the third equation, we can eliminate A and obtain 

N* = -F B 
N x = -F A 
(**) { Nf = F - AF A - BF B 
N A = F X + AF V 
^ N B — Ft + BF V . 

Conversely, the existence of a function F(t,x,<p, A, B) such that the above equa- 
tions hold clearly implies that Cn(o) — F v a G I ; but then 

C N (da) = d(C N (a)) e d(I) C I , 

and there only remains to be satisfied the condition 

4. The general isovector 
The last condition in the previous paragraph can be stated as 
(4.1) £n(P) = pa + £,da + uj/3 , 

for p a 1-form, £ a 0-form and ui a 0-form. Let D denote the coefficient of dtp in 



p ; replacing p by p — Da (which doesn't affect the validity of (4.1) as a 2 
may assume that D = 0. Setting then 



0), we 



(4.2) 

(4.3) 
and 
(4.4) 

we shall obtain a system of ten equations in F ; we shall then eliminate Ri,...,Rq. 

Identifying the coefficients of, in that order, dtdx, dtdtp, dtdA, dtdB, dxd<p, dxdA, 
dxdB, d<pdA, dipdB and dAdB yields the following system : 



p = R 1 dt + R 2 dx + R 3 dA + R 4 dB , 
t = R 5 , 



Rq 



rN v - fN A - N B + (rp -fA- B)N* + (rip -fA- B)N* - -a 2 ^ 
(4.5) = -ARi + BR 2 -(B + fA- rp)R 6 



(4.6) 



(rp - fA - B)N* - -a 2 N A = R 1 



(4.7) (rp -fA- B)N% - \a 2 N A - \a 2 N\ = BR 3 - ±* 2 R e 
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(4.8) 
(4.9) 

(4.10) 
(4.11) 



{rip -fA- B)N X B - -a 2 N A = BR t + R 5 



[B + rA- rip)Nl = R 2 



{B + fA- rip)!^ - -a 2 N l x = AR 3 + R 5 



(B + rA- rip)N B = AR± 



(4.12) 
(4.13) 



-o 2 N% = -R 3 



= -R 4 



(4.14) 



\a 2 N% 



. 



Equation (4.13) gives R4 = ; then (4.8) defines R$. Equation (4.14) is equiv- 
alent to N B = ; if that be the case, then (4.11) holds automatically. Now (4.9) 
defines R2, (4.12) defines i? 3 , (4.6) defines R\ and (4.7) defines Rq. We are left 
with equations (4.5) and (4.10) and the condition N B = 0. 

Let us begin with the last mentioned ; it is equivalent to Fbb = : F is affinc 
in B, i.e. F = c + Bd, where c and d depend only on (t, x, <p, A). Now (**) can be 
rewritten as 



(***)< 



N* = -d 

N x = -c A - Bd A 

= c - Ac A - ABd A 
N A = c x + Bd x + Ac v + ABd v 
N B =ct + Bd t + Bc v + B 2 d ip . 



From (4.12) follows 
(4.15) 

and (4.8) yields 
(4.16) 

Now we can rewrite (4.10) as 
(4.17) 



R3 



1 2 . 
--a d v , 



R 5 = (B + rA- rtp)d A - -a 2 (d x + Ad v ) . 



(B + fA-r<p)(-d A )- -<r 2 (-d x ) 



(B + rA- rip)d A - -a 2 (d x + Ad v ) . 
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Comparing the coefficients of B on both sides gives — cLa = cIa, that is cLa = 0, i.e. 
d depends only upon (t,x,ip). Now (4.17) becomes 

(4.18) \cr 2 d x = ~a 2 (d x + Ad ip ) 
that is 

(4.19) 2d x = -Ad v . 
Differentiating with respect to A leads to 

(4.20) = 2d Ax = -d v , 

whence d v — and d x = — ^Ad v = : d depends only upon t. Then (* * *) 
becomes 

( N* = -d 



(* * **) < 



N x = -c A 
N v = c - Aca 



tA 

tB 



N A = c, + Ac v 



_ N B = c t + Bdt + Bc v , 

the new unknowns being a function c(t, x, ip, A) and a function d(t), and we still 
have to satisfy equation (4.5). Now (4.9) implies R2 — 0, and (4.6) gives 

(4.21) i?i = -(rip -fA- B)c Av - ^ 2 (c xip + Ac w ) . 

Equation (4.7) now becomes 

(rip -fA- B)(-c A a) - \^ 2 (c xA + c v + Ac Aip ) - ^<J 2 (-d t ) 
= -\° 2r s > 

(4.22) 
that is : 

2 

(4.23) R e = -dt + c xA + c v + Ac Av + -^(B + fA - -rp)(-c AA ) ■ 
Eliminating R\, R 2 and i? 6 turns (4.5) into 



r(c - Ac A ) - f(c x + Ac v ) - (c t + Bd t + Bc v ) + (rip -fA- B)(-c Ax ) 
+ (rip -fA- B)(-d t ) -^a 2 (c xx + Ac vx ) 

(4.24^= -A(-(r<p -fA- B)c A , v - ^(r 2 (c xlf + Ac vv )) 

2 

-(B + fA- rip)(-d t + c xA + c v + Ac Av + —(B + fA- r<p)(-c A A)) ■ 

Both members of (4.24) are second-order polynomials in B ; equating the coeffi- 
cients of B 2 gives caa — 0, whence c is affine in A : c = e + Af with e and / 
functions of (t,x,ip). 

Equating the coefficients of B yields 

(4.25) - d t - c v + c Ax +d t = -Ac Av - (-d t + c Ax + + Ac Alfi ) 
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that is : 

(4.26) dt = 2c Ax + 2Ac Aip = 2f x + 2Af v . 

Differentiating the last equality with respect to A gives f v = (that is, / is a 
function of (t,x)), and then we get that d t = 2f x , i.e. 

(4.27) / = \d (t)x + v(t) 

for some function fi of t alone. 

Equating the constant terms in B now gives us 

re - f(e x + Af x + Ae v ) - (e t + Af t ) + {rip - fA)(-f x ) 
+ {r<p - fA){-dt) - \^ 2 {e xx + Af xx + Ae vx ) 

= A ^ 2 (e x <p + Ae w ) - (fA - rtp)(-d t + f x + e v ) . 

(4.28) 

Both sides of equation (4.28) are polynomials in A with coefficients depending 
only upon (t,x,<p). Identifying the terms in A 2 leads us to e vv = 0, whence 
e = g + hip with g and h depending only upon (t, x). The unknowns are now d and 
\x (functions of t alone) and g and h (functions of (t, x)). 

Identifying the coefficients of A on both sides of (4.29) yields 

-f(f x + e<p) - ft + ?fx + rdt - \^ 2 {j xx + e vx ) 
= ^cr 2 e xip - f(-d t + f x + e v ) , 

(4.29) 
that is 

(4.30) - f t - \o 2 f xx = o 2 h x - ff x 
or 

(4.31) a 2 h x = f^d\t)-^d"{t)x-^{t) , 
that is 

(4.32) h(t, x) = ^d (t)x - d -^x 2 - + k(t) , 

for some function k of t alone. 
The constant term gives 

(4.33) re - fe x - e t - rtpf x - ripd t - ^<J 2 e xx = rtp(-d t + f x + e v ) . 
According to the relation d t = 2f x , this becomes 

(4.34) rg - fg x ~ riph x - g t - iph t - ripd t - \^J 2 g xx - ^a 2 (ph xx = . 

Equating the constant term (in <p) of (4.34) to zero gives that g is a solution of 
(&)• 
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The value of the term in tp means that 
(4.35) - fh x - h t - rd t - \^ 2 h xx = , 

that is 

(4.36) 

This is a polynomial equation in x with functions of t as coefficients. Considering 
the coefficient of x 2 gives d (t) = 0, that is d(t) — Cit 2 + C 2 t + C3 for constants 
Ci, C 2 and C 3 . 

Equating the terms in x leads to 

(4.37) ^d(t)-^d"(t) + ^^0, 
that is 

(4.38) //'(£)= 
or 

(4.39) /Li(i) = C 4 t + C 5 

(C4, C5 real constants). 

We are left with the constant term in x : 

(4.40) - ^d (t) + f k ' (t) - r j it) + iM = , 

or 

s 2 „ . _u(i) rf'(i) ^ 



(4-41) 

for some constant C%. Therefore 



(4.42) 
and 



f(t,x) = ^d\t)x + ^(t) 

= ^x(2dt + C 2 ) + C 4 t + C 5 



h(t,x) = -I- x {2C 1 t + C 2 )~ ^x 2 

+ r t. 1 V C 6 ■ 



(4.43) 
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We have established 

Theorem 4.1. The general isovector N for (£ 2 ) is given, in terms of an arbitrary 
solution g of (£ 2 ) and six arbitrary real constants C\,...,C§, by the formulas 

N f = -dt 2 - C 2 t - C 3 
N x = - X -x(2C x t + C 2 ) - {Cit + C 5 ) 

N* = 9 + y( _L, (2Clt + C2) _^_^ 

2<7 2 (J 4 

^ = fo + ^(2C lt + C a )-^-^) 
+ |(2dt + C 2 ) 
+ A(^(2dt + C 2 )-^x 2 

. (C 4 f + C 5 ) ^ 2dt + C 2 
+ r + +C 6 ). 



N* = gt + ^-^(2C lt + C 2 ) + ^ + ^) 
+ A(xC 1 +C i ) + B(2C 1 t + C 2 ) 

_ (C 4 t + C 5 ) ^ 2dt + C 2 
+ r -5 + + C 6 ). 



5. The Lie algebra 

For u a solution of (£ 2 ), let N u denote the isovector defined by g = u and 
C\ = ... = C*6 = 0, and, for 1 < i < 6, let A 7 , denote the isovector defoned by by 
g = 0, Ci = 1 and Cj = for j 7^ i. 

The function g and ft- are determined by the isovector N : 

g = N^- <p— 
dip 

and 

h = ~a — ; 

dip 

we shall denote them repectively by g^ and fiN- As seen in §4, is a function of 
(t,x). 
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Lemma 5.1. For all (M, N) G Q 

- M * d9N -X- M x d9N -L n h AT* 9 9M , 
9[M,N] — M -gj- + M -q^- + 9MIIN - N JV gNflM 

and 

h l\jt dh N , n/r x dh N AT t dh M ATX 0h M 

h [M , N] -M—+M —-N—-N —. 

Proof. One has 

[M, N] v = M(N V ) — N{M V ) 

= M(g N + fh N ) - N(g M + iph M ) 

ot ox at ox 

+M*h N -(N^ + N^ 
ot ox 

+^N^ + N^)+N*h M ) 

ot ox ot ox 

+ (5m + iph M )h N - N — N ~q~T 

~ dt~ ~dx~^ _ ^ 9N lfhN ^ hM 

- I m * dgN j. m x dgN m n h yt d 9M vx dg M , . 

the result follows. □ 
Let us set 



and 



n = {N G G\g N = 0} 



J = {N G Q\hpf = 0} 



Proposition 5.2. J7" is an idea/ o/ Q and % is a subalgebra of Q . Furthermore 
Q = H (B J, FL has dimension 6 and admits (N\, Nq) as a basis, and J = 
{N u \u solution of (£2)}- 

Proof. Clearly, N G J if and only if C\ = ... = C% = ; in particular, if TV G J then 
TV* = N x = 0. Therefore M G Q and iV G J imply /i [MjjV ] = 0, i.e. [M, N] e J: J 
is an ideal of Q. 

Furthermore, gu = gN = imply g[M,N] = : % is a subalgebra of <J. The last 
two assertions clearly hold. □ 

6. Some symmetries 

Let N G G, let k G R, and let ^ be a solution of (£ 2 ) ; then e nN maps 
(i, x, <p, A, B) to (t K ,x K ,(p K , A K ,B K ) ; setting 

(6.1) p K = 4> K {t K ,x K ) , 
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it follows that is also a solution of {£2)- Wc shall denote 
(6.2) e K * : <p ^ ^ 

the associated one-parameter group. 
Lemma 6.1. If 



N = N 1 — + N x — + NV— +...66 
ot as atp 



then 



that is 



N(<p) = (dt 2 + C 2 t + C 3 )^ + (\x(2C\t + C 2 ) + (C 4 t + C 5 ))^ 

+g + v (^ x (2C 1 t + C 2 ) - ^x 2 - ^ - ^(Cit 2 + Czt + C 3 ) 
, „C 4 t + C 5 2dt + C 2 

+r a 2 + 4 + • 

Proof. Let us rewrite (6.1) as 

(e KN (<p))(t,x) = (e^(^)(e« N (t),e* N (x)) . 
Developping at order one in k gives 

<p + nN v = (p(t + nN\x + KN x ) + KN(Lp)+o(K) 

= y + kN*^ + kN*^ + ltN((p) + (k) , 

ot ox 

whence the result. □ 
We shall set 

(6.3) C^(t,S):=(e^ V )(t,ln(S)) 
and 

(6.4) (•■: : 
(1 < 3 < 6). 

Some of these transforms can actually be explicitly computed ; for instance 

(6-5) 
whence 

~2 

(6.6) (e Kl * 3 <p)(t,x) = e~^(p(t + K,x) 
and 

(6.7) C 3 K (i,5) = e^C(t + K ,S). 



(6.8) 
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therefore 

Ki . _ . Kj t 

(6.9) (e K %)(t, x) = e^ (r< ~ ~ 2o^ ' <p{t, x + nt) 
and 

(6.10) C£(t,S) = e2^ ( r_K ^ _ ^C(i,e Kt ^) . 



whence 



(Ml) 



(6.12) (e KN5 ip)(t,x) = ev 2 ip{t,x + k) , 
and 

(6.13) C 5 K (i,S) = e^C(t,e K S*) • 
Lastly, 

(6.14) N 6 (<p) = <p 
whence 

(6.15) (e Kl ^ e ip)(t,x) = e K ip(t,x) 
and 

(6.16) Cg(t,S) = e K C(t,S) . 

Corollary 6.2. Let C denote a solution of (£). Then, for each ceR, the (3 < 
i < 6) defined by (6.7), (6.10), (6.13) and (6.16) are also solutions of (£). 

Together the transformations described in (6.7) and (6.16) come from the in- 
variance of the original equation via multiplication of the solution by a scalar and 
translation in time ; (6.13) then comes from the homogeneity in S — these could 
be expected. The transformation given by (6.10) is, however, not so easy to under- 
stand. It would be interesting to find some financial interpretation for it. 
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